Selective coupling of superconducting qubits via tunable stripline cavity. 
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We theoretically investigate selective coupling of superconducting charge qubits mediated by 
a superconducting stripline cavity with a tunable resonance frequency. The frequency control is 
provided by a flux biased dc-SQUID attached to the cavity. Selective entanglement of the qubit 
states is achieved by sweeping the cavity frequency through the qubit-cavity resonances. The circuit 
is able to accommodate several qubits and allows to keep the qubits at their optimal points with 
respect to decoherence during the whole operation. We derive an effective quantum Hamiltonian 
for the basic, two- qubit- cavity system, and analyze appropriate circuit parameters. We present 
a protocol for performing Bell inequality measurements, and discuss a composite pulse sequence 
generating a universal control-phase gate. 



o 
o 

u • 

Oh 



X3 
O 

> 

O 
(N 
00 

o 

O 

-I— > 



I 

O 

o 



X 



Coherent coupling of superconducting qubits has been 
experimentally demonstrated for all major qubit types 
(chargei^, flus^*^, and phased qubits) using perma- 
nent direct qubit-qubit coupling, capacitive or inductive. 
A major challenge is to implement a tunable coupling of 
qubits required for any useful gate operation. Numerous 
suggestions in this direction have been discussed in recent 
literature together with related quantum gate protocols 
(for a review see, e.g. Ref. 

There are two conceptually different approaches to the 
tunable coupling. The first approach is to employ di- 
rect coupling schemes using Josephson junctions in the 
non-resonant regime as passive controllable elements, ei- 
ther capacitive,^ or inductiveiiSiiiiiSiiiii The second ap- 
proach, which we adopt in this paper, suggests qubit 
coupling via a dynamic intermediate element, e.g., LC- 
oscillator or Josephson junction, which becomes entan- 
gled with a qubit during a two-qubit operation. In this 
scheme, the entanglement is achieved by tuning the qubit 
and the mediator in resonance, and then transferring the 
entanglement to another qubit by tuning the mediator 
and the second qubit in the resonance. Such coupling 
method has been first suggestedi^ and experimentally 
testediS for the ion trap qubits. For superconducting 
qubits, qubit-oscillator entanglement has been demon- 
strated experimentally for a charge qubit coupled to a 
microwave stripline cavity^ii and a flux qubit coupled to 
a SQUID oscillator fiSJA the gate protocols based on con- 
trollable qubit-osc illator coupling have been theoretically 
discussed in Refs. l20l2lL 

The experimental setup with the qubit coupling to a 
distributed oscillator - stripline cavitjii2i2^ possesses po- 
tential for scalability - several qubits can be coupled to 
the cavity. In this paper we investigate the possibility to 
use this setup for implementation of tunable qubit-qubit 
coupling and simple gate operations. Tunable qubit- 
cavity coupling is achieved by varying the cavity fre- 
quency by controlling magnetic flux through a dc-SQUID 
attached to the cavity (see Fig.0. An advantage of this 
method is the possibility to keep the qubits at the opti- 
mal points with respect to decoherence during the whole 
two-qubit operation. The qubits coupled to the cavity 
must have different frequencies, and the cavity in the 
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FIG. 1: Sketch of the device: charge qubits (single Cooper 
pair boxes, SCB) coupled capacitively (Cc) to a stripline cav- 
ity integrated with a dc-SQUID formed by two large Joseph- 
son junctions (JJ); cavity eigenfrequency is controlled by mag- 
netic flux $ through the SQUID. 



idle regime must be tuned away from resonance with all 
of the qubits. Selective addressing of a particular qubit 
is achieved by relatively slow passage through the res- 
onance of a selected qubit, while other resonances are 
rapidly passed. The speed of the active resonant pas- 
sage should be comparable to the qubit-cavity coupling 
frequency while the rapid passages should be fast on this 
scale, but slow on the scale of the cavity eigenfrequency in 
order to avoid cavity excitation. This strategy requires 
narrow width of the qubit-cavity resonances compared 
to the differences in the qubit frequencies, determined 
by the available interval of the cavity frequency divided 
by the number of attached qubits. This consideration 
simultaneously imposes a limit on the maximum num- 
ber of employed qubits. Denoting the difference in the 
qubit energies, AEj, the coupling energy, k, the max- 
imum variation of the cavity frequency, Aw^, and the 
number of qubits, N, we summarize the above argu- 
ments with relations, k <C Ai?j, TV ~ hAujk/ AEj. In 
the off-resonance state, the qubit-qubit coupling strength 
is smaller than the on-resonance coupling by the ratio, 
K/{hujk-Ej) < 1. 

In the first part of the present paper, we analyze the 
quantum electrical circuit consisting of superconducting 
stripline cavity, dc-SQUID, and single Cooper pair box 
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(SCB) qubits, and derive an effective quantum Hamil- 
tonian for this circuit, and discuss the relevant circuit 
parameters. 

Then, on the basis of the derived Hamiltonian, we dis- 
cuss the Bell measurement protocol and a protocol for 
a conditional phase gate. We consider creating maxi- 
mally entangled two-qubit states (Bell-states) by sequen- 
tially sweeping the cavity through the resonances with 
the two qubits^S and discuss the protocol for measuring 
the CHSH correlation function^'^ for such states, which 
is equivalent to testing the Bell inequality. 

While considering the universal two-qubit gate, we 
take into account an important feature of our system - the 
linearity of the cavity, which does not allow implementa- 
tion of the V SWAP gate.^^ We argue that the control- 
phase gate (CPHASE) is a genuine two-qubit gate for 
our system (cf. Ref. l2(ll) . We consider a protocol for 
this gate, which is much faster than the one suggested 
in Ref. .2(k the present one being based on the resonant 
rather than dispersive qubit-oscillator coupling. A major 
difficulty for constructing such a protocol is the genera- 
tion of the single- and two-photon states in the cavity 
(for the cavity initialized in the ground state); elimina- 
tion of these auxiliary photon states requires a complex 
pulse sequencei^ 

In this paper, we explicitly discuss the coupling of 
charge qubits; however, the method of derivation of the 
effective quantum Hamiltonian also applies, with minor 
modifications, to the flux qubits, and the quantum pro- 
tocols studied can be extended to this type of qubit sys- 
tems. 



I. CAVITY WITH VARIABLE FREQUENCY 

The resonant frequency of a ID stripline cavity de- 
pends on the boundary conditions. For example, if one 
end of the cavity is open while the other is connected 
to the ground, the spatial distribution of the supercon- 
ducting phase along the cavity has a maximum at the 
open end and a node at the grounded end. This cor- 
responds to the quarter- wavelength resonator, d = A/4, 
with the eigenmode wave vectors, fc„ = {n/d){n + 1/2), 
where d is the cavity length, the eigenmode frequencies 
being a;„ — {nv/d){n + 1/2), where v is the velocity of 
the electromagnetic waves in the cavity. If, on the other 
hand, the second end is disconnected from the ground, 
the eigenmode wave vectors become, fc„ = (7r/d)n, giv- 
ing the frequencies, uJn — {■nvld)n. The role of the dc- 
SQUID attached to the stripline cavity in Fig. ^ is to 
vary the boundary condition at the right end: the first 
case (node) corresponds to a very large (formally infinite) 
Josephson energy of the SQUID, while the second case 
(antinode) corresponds to the fully suppressed Joseph- 
son energy. Thus, ideally, by changing the biasing mag- 
netic fiux through the SQUID by half a flux quantum, 
< $ < $o/2(mod$o), (^0 = h/2e), one should be 
able to sweep the eigenmode frequencies within the inter- 
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FIG. 2: Equivalent circuit for the device in Fig. Q chain of 
LC-oscillators represents the stripline cavity, and cj>N are 
superconducting phase values at the ends of the cavity, ipj 
and are local phase values where the qubits are attached; 
attached dc-SQUID has effective flux-dependent Josephson 
energy, Ejs{f), and capacitance Cs, control line for tuning 
the SQUID is shown at the right; SCB qubits are coupled to 
the cavity via small capacitances, Cci and Cc2- 



vals, {'Kv/d)n < w„ < {-Kv/d){n -f 1/2). In practice these 
intervals are more narrow due to a finite maximum and 
non-zero minimum Josephson energies of the SQUID. 

For a given eigenmode, the integrated stripline -I- 
SQUID system behaves as a lumped oscillator with vari- 
able frequency. Our goal in this section will be to derive 
an effective classical Lagrangian^^'^^ for this oscillator. 
To this end we consider in Fig.|5]an equivalent circuit for 
the device depicted in Fig. ^ A discrete chain of identi- 
cal LC-oscillators, with phases across the chain capac- 
itors (i=l,. . . ,N), represents the stripline cavity; the dc 
SQUID is directly attached at the right end of the chain, 
while the superconducting Cooper pair boxes (SCB) are 
attached via small coupling capacitors, Cci and Cc2 to 
the chain nodes with local phases, <j)j and (pi (for simplic- 
ity we consider only two attached SCBs). The classical 
Lagrangian for this circuit, 

— ^SL ~t- ^squid ^" ^ ^ ["^Q^i ^ ^couplj] ; (1) 
J = l,2 



consists of the stripline Lagrangian, 
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the SQUID Lagrangian, 
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the Lagrangians of the SCBs, 
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and the capacitive SCB-stripline coupling, 

^2 {4>J + <Pc 



2(2e)2 



(5) 



The SQUID junction variables are related through the 
flux quantization relation, (j)s.i — 0s, 2 = /, to an exter- 
nally applied magnetic flux, $ — (<i>o/27r)/ = {h/2e)f, 
threading the SQUID ring. The self inductance of the 
SQUID ring is assumed to be negligibly small compared 
to the Josephson inductances of the SQUID junctions. 
Then the SQUID can be described as a single junc- 
tion with effective capacitance, Cs, and flux-dependent 
Josephson energy. 
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A. Linear approximation 




FIG. 3: Solution of dispersion equation IIH for first mode, 
kd < 77 /2 {d < A/4), for large (a) and small (b) Josephson 
energies of the SQUID {{2e/hf L^avEs = 16 and 4, respec- 
tively); inset shows corresponding spatial distributions of the 
phase (f)/(j>{0) in the cavity. 



Let us assume small amplitude of the plasma oscilla- 
tion in the SQUID, 0s ^ 1, which implies the phase 
regime for the SQUID, Ej^if) > (2e)V2Cs, and then 
adopt the harmonic oscillator approximation in Eq. Q , 
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where (ps — {4>s,i+4's,2)/'2+'r]{f), r]{f) is a constant phase 
shift, which can be neglected for adiabatic flux varia- 
tions. The SQUID Josephson energy Ejs{f), Eq. 10, 
reaches its maximum at zero magnetic flux, EJ"^^ = 
Ejs,i+Ejs^2i while the minimum is approached at / = tt: 
^rnjn ^ \Ej,^i-Ejs^2\ with E^™™ > duc to the SQUID 
asymmetry. 

To proceed to a continuum description of the stripline 
cavity, we introduce the distance Aa; between nodes i 
and i + I, and express the stripline Lagrangian, Eq. |(2Jl, 
in terms of the stripline capacitance and inductance per 
unit length, 

Co = C/Ax, Lq = L/Ax. 

Let Ax go to zero and transform the node index i into 
the continuous variable x. In the bulk of the cavity, the 
equation of motion of the field is a wave equation, 



{x,t) -v'^(j)"{x,t) = 0, 



(8) 



where v = I/^/LqCq is the wave velocity. It is con- 
venient to express the wave velocity through the cav- 
ity inductance, Lcav = dL^, and the cavity capacitance, 
Ccav = dC^, 



The boundary condition at the cavity open end (x 



0), 



0'(O,t) = O, 



(9) 



requires that the only allowed solutions are of the form 
(\)(x, t) — <pi sin(fcwt) cos(fcx). The boundary condition at 
the cavity right end (x = d) reads, 
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^'{d,t)+Ejs{f)cP{d,t)^0. 



A dispersion equation for the cavity eigenmodes results 
from Eq. (|l()|l using the bulk solution to Eq. JHJ, and 
takes the form, 



(kd) tan(fcd) = ^ L,avEs (/) 



a. 



(kdr. (11) 



The solutions to this dispersion equation form an infi- 
nite set of eigenmodes with wavelengths A = 27r/fc and 
frequencies uJk — kv. 

The solutions to Eq. Hll|l are illustrated in Fig. O 
they are given by the intersection points of the function, 
kdtSLii{kd), with the parabola, which is almost flat in the 
practically relevant limit, Cs/Ccav ^ 1- The zeros of this 
function {kd = wk, = 0, 1, ...) correspond to an open 
right end of the cavity (disconnected SQUID), while the 
singular points (kd — it /2 + mi) correspond to a closed 
cavity end (short circuited SQUID). These limits of the 
variation of the cavity wave eigenvectors, mr < knd < 
7r/2-|-7rn, can be achieved when Es{f) varies between oo 
and 0; thus ideally the frequency can be tuned between 
mt/y/LcavCcav and (7r/2 + Tin) / LcavCcav In practice, 
the available frequency range is smaller, being limited by 
the value of the parameter [(2e)^/ {2Ti^)]LcavEs{^), which 
should be chosen large, and the minimum value of the 
SQUID Josephson energy. Eg (tt) , allowed by the SQUID 
asymmetry. 

Let us return to the Lagrangian of the stripline cavity 
and the SQUID, Eqs. l(2Jl and ||7J), and only consider a 
single eigenmode, (j){x,t) = (f)i{t) cos(kx). In the contin- 
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uum limit the Lagrangian will then take the form, 



dx 



+ 



h 

2e 
2(2e)- 



Co4>lcos'^{kx) 0?fc2sin^(fca;) 

2 QjTn 



After performing integration over x, and using the disper- 
sion equation we arrive at the effective LC-oscillator 
Lagrangian representing the integrated cavity+SQUID 
system, 



r 

*^os 



n^Ck 

2(2e)' 



2(2e)2Lfc 



^l{t). (13) 



The oscillator is described by the effective fc-dependent 
capacitance Ck, 



Ck 



Cc. 



^ ^ sin(2fcd)\ 



2 V ' 2fcd 
and the effective inductance L^, 



+ CsCOs'^{kd), (14) 
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Lk 



{kdf 
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sin(2fcd) 2Cs 



2kd 



cos^(fcd) .(15) 



The frequency of the effective oscillator, defined in the 
usual way, ujk = \/ ^/CkLk^ is equal to the frequency of 
the chosen cavity eigenmode as one should expect. 



kd 



i^k 



^ ^cav Ccav 



(16) 



B. Non-linear correction 

While the stripline cavity alone is a linear electromag- 
netic system, attaching the dc-SQUID makes the inte- 
grated system non-linear. Non-linearity will introduce a 
non-equidistant correction to the quantized energy spec- 
trum of the cavity, which may affect the gate protocols; 
in particular it is harmful for the conditional phase gate 
protocol considered later in the paper. Therefore it is 
important to estimate non-linear effects produced by the 
SQUID on the cavity. 

To this end we expand the SQUID potential in the 
boundary condition in Eq. I|10|l . assuming 0^ <C 1, and 
keep a small cubic term, 



(2e)2 



kd,t) 



h^d 



(2e)2L, 



■(f>'{d,t) 



(t>{d,t) - - ^^d,t) 
6 



0. 



(17) 



The cubic term will introduce the third harmonic in the 
cavity. 



(x, t) = 01 sm{kvt) cos{kx) + 4>z smi^kvt) cos{ikx) 



(18) 



whose amplitude 
condition, H17|l. 



can be found from the boundary 



»3 = 



Ak 
24 



(19) 



^lcos'^{kd)-^4p-'t>lcos^{kd). (12) Ak 



Lul cos^ {kd) 



[^2 _ 9{kv)^] cos(3fcd) - {ikd/L^avCs) sin(3fcd) ' 

Here we introduced the plasma frequency of the SQUID, 

ujl{f)^{2e/hf[Ej,{J)/CA. (20) 

The cubic term also produces a shift of the resonance 
frequency given by the corrected dispersion equation. 



[kd) tan(fc(i) 



{2ef 



Cca 



Esif) 
{kdf . 



1 



1~ -ujj. cos {kd)4> 



(21) 



Taking the relation (|11|) into account, and omitting a 
small term ~ Cg/Ck ^ 1, we obtain the relative shift of 
the frequency. 



5k 
T 



Bk4>l- 



Bk 



(1/4) cos^ {kd) 
1 + 2kd/ sin(2fcd) 



(22) 



Such an amplitude-dependent frequency shift, on the 
other hand, can be recovered from the effective oscilla- 
tor Lagrangian in Eq. (|13|l by adding the following non- 
quadratic term, 



2{2e)^Lk 



Bk<j>\. 



(23) 



In the quantum regime, such a term will produce a devi- 
ation from the equidistant energy spectrum of the cavity. 
The magnitude of this deviation in the first perturbative 
order reads, 



= - 



Bk 



2{2e)^Lk 
+ 6n + 3 



■ BkEck, 



(24) 



where Eck — {2ef' /2Ck is the charging energy of the 
cavity, and n is the energy level number. Thus we see 
that the non-linear effect is proportional to the charging 
energy of the cavity. In order to neglect the non-linear 
effect, this energy must be much smaller than the energy 
of the qubit coupling to the cavity (see below). 



II. QUBIT COUPLING TO THE CAVITY 

Now we take the SCB qubits in Eq. into the 

consideration, assuming that the coupling capacitances, 
Ccj are small enough that the perturbation of the cavity 
eigenmodes due to the SCBs is negligible. 
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The cavity field at the point where SCB j is coupled 
is related to the effective oscillator variable (pi by the 
relation, = (f)iCos{kxj), where Xj is the position of 
qubit j along the cavity. The coupling is described by 
the cross term in Eq. Q, 



int.j 



(2e) 



■ Cc,j C0s(kXj)(j)l4>q,. 



(25) 



the quadratic terms in Eq. Q give small renormalization 
of the qubit capacitance, Csj — Cj + Cc.j + Cg, and the 
oscillator capacitance. 

We transform the capacitive interaction into an induc- 
tive form, 

Ant,i = ctjEjjipi sin ipqj, (26) 
using the transformations^, 

Cj)q,j ^ (l)q,j + a^jCpl, (27) 

with the coupling constant, 

aj = Cc.j cos{kxj)/Csj, aj ^ 1. (28) 

The SCB Lagrangian does not change during the trans- 
formation, whereas the oscillator undergoes displace- 
ment, 



^LC — 



h Ck -,2 



i'l - ■;^C'g(Q;iVg,i -f a2Vg^2)4>i 



(29) 



and small renormalization of the effective capacitance, 
Ck -> Cfc - Y,j=i,2 ^c.j co8{kxj)[cos{kxj) - aj]. 



Effective Hamiltonian 



We obtain the classical circuit Hamiltonian, 



i=i,2 



nt.,j\ ~t~ Hqsci 



(30) 



by introducing the conjugate momenta n = {l/fi)dL/d(f>i 
and rij — {l/h)dL/d(t)qj. Each SCB is described by the 
Hamiltonian, 



Hj — EcjijLj — Ugj)"^ — Ejj cos( 



9J' 



(31) 



where Ecj — (2e)^/ (2Cs j) is the charging energy of the 

CgVgj / {2e) is the (dimension- 



SCB island, and 

less) charge on the island induced by the gate-voltage. 

According to Eq. H29|l . the gate voltages also induce 
a charge, uq — aiUg^i + a2"-g,2, on the oscillator. Be- 
cause the oscillator charge is not quantized, this induced 
charge does not have any physical meaning, and can 
be eliminated using the gauge transformation WHoscU 



with U = exp(— mo0i). The oscillator Hamiltonian then 
readsS^, 



Hose = EcklT? + ELk<Pi, 



(32) 



where Eck — (2e)^/(2Cfe) is the charging energy of the 
oscillator and E^k — ?i^/(2(2e)^Lfe) its effective inductive 
energy. The interaction term in Eq. (|30|l is given by the 
expression in Eq. (|26|) with the opposite sign, H int.j — 

The Hamiltonian H3U|I is quantized by imposing the 
canonical commutation relations, [(f)qj,nk] = iSjk, 
[4)i,n\ = i. For later convenience, the oscillator is de- 
scribed in terms of the ladder operators. 



with [a, a^] = 1. The quantized oscillator Hamiltonian 
then reads, 



Hose = TT-i^k i^ + a^a 



(33) 



with LOk given by Eq. (|16() . 

The Coulomb blockade effect in the SCB is taken into 
account by considering the periodicity of the SCB po- 
tential and imposing 27r-periodic boundary conditions on 
the wave function with respect to the phase (j)qj- The re- 
sult is charge quantization on the island. Assuming the 
charge regime, Ec ^ Ej for the SCB, and keeping the 
system at low temperature {ksT <^ Ec) and close to the 
charge degeneracy point Ug = 1/2, restricts the number 
of excess charges on the island to zero or one Cooper 
pair. This allows us to truncate the SCB Hilbert space 
to these two lowest charge states |0) and |1) . 

It is advantageous to operate at the qubit charge 
degeneracy point, where the decoherence effect is 
minimized , and allow only small departures, Sugj ~ 
1/2 — ngj", from this point during single qubit operations. 
Considering this, we write the quantized qubit Hamilto- 
nian in the qubit eigenbasis at the charge degeneracy 
point, [|g),|e)] = [|0) + |l),|0)-|l)]. 



H] = E, 



c,jOng^jax^j 



(34) 



The interaction, Eq. H26(l . is proportional to sini^qj , 
which transforms into ay.j during the quantization proce- 
dure. It will be helpful during the discussion of two-qubit 
operations to express the interaction through the raising 
(lowering) qubit operators <J-\e) = \g), a+\g) = |e), 

Thus the quantized interaction Hamiltonian reads, 



H 



int,j 



i^(a + at)(f7_,, -a+,j), (35) 
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where the interaction energy Kj is determined by the cou- 
pHng constant Uj in Eq. (|28|l . 

''^^'^'^'''[j§t) ■ ^^^^ 

Equations H35II . H3(j|l were derived for the charge hmit, 
Ec 'S> Ej. However, they remain vahd quahtatively also 
in the charge-phase regime, Ec ^ Ej, which is more 
advantageous from the point of view of decoherence, as 
is weU established.'^^ In this regime, the lowest Bloch 
states of the SCB Hamiltonian rather than the charge 
states form the computational basis. This is fully consis- 
tent with the quantum capacitance readout method'^^ for 
our system, which realizes projective measurement on the 
qubit eigenbasis. Transformation of the SCB eigenbasis 
from the charge regime to the charge-phase regime with 
increasing ratio Ej / Eq was analyzed in Ref. l33i Apply- 
ing this analysis to the present case in the Appendix we 
find that the qubit-cavity coupling remains transversal 
acquiring the form in Eq. ljA12|l . 

Hrnt.j = — ^ — {'^]^+,J + 1 (37) 

with the coupling constant kj differing from the coupling 
Kj of the charge regime in Eq. (|36|l by a complex numer- 
ical function f{Ejj/Ec,j) ^ 1, kj = f{Ej^j/ Ec,j)K.y 
Such a modification does not change qualitatively the 
resonant qubit-cavity dynamics discussed in the next sec- 
tions. 



B. Constraints 

To conclude our discussion of the qubit-cavity circuit, 
we summarize the imposed constraints on the circuit pa- 
rameters required for a proper functioning of the circuit. 

First, we required the phase regime for the SQUID 
and the cavity, implying (2e)^/(2Cs) — Ecs <C hujs <C 
Ejs{f), and Eck ^ ^^fc <C E^k, respectively. The cav- 
ity capacitance in practice greatly exceeds, by several or- 
ders, the SQUID junction capacitances, Ck ^ Cg, while 
the cavity inductance must be comparable to the SQUID 
variable inductance, which is required by the dispersion 
equation Hll() . for kd ^ 1. Thus, the SQUID plasma fre- 
quency is typically much larger than the cavity frequency, 
LOs S> cjfe. The latter, in its turn, must be comparable to 
the qubit frequencies to provide the resonant coupling, 
Wfc ~ Eji/fi,Ej2/h. 

The qubit interaction with the cavity must not be too 
strong in order to provide sufficient off-resonance decou- 
pling, K - {Cc/Cs)Ej < \Ej^2 - Ej^i\. On the other 
hand, it must exceed the variation of the level spacings 
in the cavity energy spectrum, Eq. (|^ . k Eck, caused 
by non-linearity. 

All these requirements can be collected in a chain of 
inequalities formulating the hierarchy of relevant circuit 



energies, 

Eck<^n < \Eji ~ Ej2\ < Tiujk ^ Ej, 

hivk, Ecs < huJs <^ ELk Ejs. (38) 

There is an additional requirement imposed on the 
lower bound of the variation of the Josephson energy 
of the SQUID: the critical current through the SQUID 
should be much larger than the amplitude of the current 
fluctuations in the cavity. The latter is estimated for the 
zero point fluctuations as {h/2e){kd/ Lcav){Eck/ E^kY^^^ ^ 
while the critical current of the SQUID is {2e/K)Ejs- 
Thus the requirement is equivalent, by virtue of Eq. Hll|l . 
to the inequality tan(fcd) > {Eck/ E^kY^^ ■ 

III. TWO-QUBIT OPERATIONS 

A general way of performing two-qubit operations is 
to sequentially drive the cavity frequency in and out of 
resonance with the qubits, i.e. TiiOk — Ejj = S, \S\ ^ Kj. 
The speed should be comparable to the scale of the qubit- 
cavity coupling frequency but small on the scale of the 
oscillator frequency, to prevent unwanted excitation of 
higher oscillator levels, dtOJkl^k ^ i^/^ ^ hujk- Dur- 
ing the two-qubit operations, each qubit is parked at its 
charge degeneracy point, and an appropriate difference in 
the qubit energies, \Ej^i — -Ej,2| ^ k, prevents the oscil- 
lator from simultaneously interacting with both qubits. 

Consider, as an example, qubit 1 in resonance with 
the oscillator. Discarding fast terms (^g^±^i'^k+Ej^i/K)t ^^^^ 
^±i{ujk±E.,,2/n)t-j interaction picture (with Hq = 

Hi + H2 + Hose), which average to zero on the time scale 
of the qubit-oscillator interaction in the rotating wave ap- 
proximation (RWA), the qubit-oscillator interaction term 
reads. 

Hint « Y ^ • (39) 

The only non-zero interaction matrix elements, in the 
qubit 1 - oscillator basis, are (e, nja+alg, n + 1) = (g, n + 
l|fT_a^|e, n) = ^/n + 1 between the levels which are close 
to being degenerate, En+i,g — En^e = S. The (truncated) 
Hamiltonian of this subspace reads, 

H = - ^a, + '^V^^cTy. (40) 

While discussing the entangling operations, in the fol- 
lowing we assume, to obtain analytical results, rectangu- 
lar pulse shapes bringing the cavity in and out of exact 
resonance with the qubit (^ = 0). The result of such an 
operation is given by the unitary matrix, 

= ( ^ose - sine \ ^ ^ 

^ ' \ smt' cost' / ' 2 ' \ ^ 

where T is the pulse duration. A slight detuning within 
the allowed interval 5 <C Kj and smoother pulse shape 
will not qualitatively alter the protocols. 
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A. Bell measurement protocol 

Maximally entangled two-qubit states (Bell-states) are 
constructed from an initial state with one excitation, say 
\egO), and partly moving this excitation to the other 
qubit by sequentially sweeping the oscillator through res- 
onance with both qubits.^° The oscillator starts and ends 
in the ground state. The first pulse brings the cavity to 





U2 CI' 2) 




M 












Ui('1>i) 


/A 















T, 



FIG. 4: Protocol for creating a Bell-pair: the cavity frequency 
is sequentially swept through resonances with both qubits; at 
the first resonance the oscillator is entangled with qubit f , at 
the next resonance the oscillator swaps its state onto qubit 
2 and ends up in the ground state. A Bell measurement is 
performed by applying Rabi pulses to non-interacting qubits, 
and projecting on the qubit eigenbasis, {\g), |e)}, by measur- 
ing quantum capacitance. 

resonance with the excited qubit 1, hujk = Ei (see Fig.^ 
during a time satisfying the relation, kiTi ~ hTT/2 (7r/2- 
pulse). Then the cavity, which has become entangled 
with qubit 1, 



1 

71 



{\egQ) + \ggl))., 



is driven towards resonance with qubit 2. The accumu- 
lated phase during the free evolution {Tiujk 7^ Ei,E2 for 
a time T2) is ip2, 



1 



-=(|egO)-t-e'^^|55l», ^2 = 

Jo 



^1 u\ 



dt. 
^(42) 

After the oscillator has reached resonance with qubit 2, 
huJk = E2, it stays in resonance during the time, K2T3 = 
Tltt (TT-pulse). The system evolves to the state. 



V2 



Choosing the time T2 of free evolution such that the ac- 
cumulated phase, (p2, satisfies the equation, exp[i(£'i — 
E2)T3 / h + iip2] = Til the two-qubit system ends up in 
one of the Bell states, 



|vl/±) = i=(|e5)±|ffe))|0). 



(43) 



Note that the cavity has returned to its ground state. 



To perform a test of the non-classical statistical prop- 
erties of the Bell state, Eq. (|43|l . it is convenient to im- 
plement a protocol similar to the one of Ref. |^ for mea- 
suring the CHSH inequality^''. The protocol consists of 
two independent rotations of the uncoupled qubits by 
means of applying Rabi pulses, and then measuring the 
quantum capacitances of both SCB^i. The latter realizes 
a projective measurement on the qubit eigenstates, \gj), 

The resonant 7r/2-pulse applied to the gate of the j-th 
qubit. 



^^gjit) = cos{Ejj t), 



(44) 



during the time Ta such that Ta = hn/ {2EcSn'^ j) pro- 
duces a unitary transformation. 



1 



1 



■',a = EjjTa/h, 



(45) 

By applying such pulses to both qubits, one gets the 
state, 



*±('/'l>2) = C/i((/.i)[/2(02)* 



(46) 



After the measurement, the qubits will be found either in 
similar states (both qubits in the ground, \gg), or excited, 
|ee), states), or in different states {\ge) or \eg)). The 
corresponding correlation function reads, 

9±('/'l,02) = {^±i<Pl,(l)2)Wz,l(7zM'^±{(pl,(t)2))- (47) 

By repeating the procedure with a 7r/2-pulse of slightly 
different duration, T;,, with Ti, — Ta ^ h/Ejj <C Ta, we 
obtain the four correlation functions, 

qi4>la,4>2a), g(</il6, 02a), qi't'la , <j>2b) , q{4>lb , 4'2b) ■ (48) 

According to the analysis in Ref. 23, the quantity, 

B = |g(01a, (/'2a)-'?((/>16, 02a) | + k('/>la, 026) +(7(016, 02b)|, 

(49) 

has an upper bound for classical statistics, B < 2. For 
the state in Eq. (|43|l . however, we have. 



^±(01,02) = ±COS(01 - 02), 



(50) 



and the upper bound for B becomes B < 2\/2, the equal- 
ity being achieved for, 



37r 



'Jla = -92a 



•Jib = -(P2b 



(51) 



B. Two-qubit control-phase gate 



In this section we modify the Bell state construc- 
tion to implementing a control-phase (CPHASE) two- 
qubit gate. This gate has the diagonal form: |a/30) 
exp(i0Q;3)|Q;/3O) (000 = 001 = 0io = 0, 0ii = tt), and it 
is equivalent to the CNOT gate (up to local rotations). 
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To generate such a diagonal gate, we adopt the following 
strategy: first tune the oscillator through resonance with 
both qubits performing 7r-pulse swaps in every step, and 
then reverse the sequence, as shown in figure |S1 With an 
even number of swaps at every level, clearly the resulting 
gate will be diagonal. 



Ti T2 T3 T4 T5 t 

FIG. 5: Pulse sequence producing (trivial) diagonal gate; dur- 
ing time Ti, qubit 1 swaps its state onto the oscillator, then 
the oscillator interacts with qubit 2 before swapping its state 
back onto qubit 1; free evolution during time Ts is added to 
annihilate two-photon state in the cavity. 

Such a strategy would indeed produce the CPHASE 
gate provided the qubit-cavity coupling-constant in 
Eq. (|4()|l does not depend on the number of photons. In- 
deed, applying the rectangular 7r-pulse to pass the first 
resonance, KiTi = ttH, and then the 7r-pulse for the sec- 
ond resonance, K2T2 = nh, and then reversing the pulse 
sequence, K2T4 = ttH, and kiT^ = nh, see Fig. |31 we 
induce the following transitions. 



(52) 



which generate the CPHASE gate. The problem is, how- 
ever, that it is not possible to make the swaps for all 
the states at the second resonance, e.g., l^eO) Isffl), 
and \gel) —^ 1^52), with the same pulse because of dif- 
ferent values of the coupling constants, K2, and K2\/2, 
respectively. Thus the 7r-pulse for the first transition will 
necessarily produce a state superposition for the second 
transition and vice versa. There is a possibility to anni- 
hilate this superposition by inserting an interference loop 
in the pulse sequence, see Fig. ]S[ namely, by departing 
from the resonance for a while to accumulate an appro- 
priate phase shift during the free evolution; the required 
time, Ta for such an excursion is given by equation. 





IggO) ^ 


\ggO) ^ 


I.95O) ^ 


+ I.95O) 


\geO) ^ 


\geO) ^ 


l.9.9l> ^ 


-\geQ) ^ 


-\geO) 


|e<?0) ^ 




-\geO) ^ 




+ \egO) 


|eeO) ^ 


\gei) ^ 


1.9.92) ^ 


-\gei) ^ 


+ |eeO), 



73 = 



— ^k[t) 

n 



dt = TT (mod 2tt). 



(53) 



However, it is easy to check that such a protocol will 
generate a non-entangling, trivial gate. 

It turns out that by adding two more swap segments 
with interference loops to the previous protocol it is 



possible to obtain a sufficient amount of free parame- 
ters to annihilate the state superpositions, and to obtain 
an entangling gate24. The pulse sequence is shown in 
Fig. IHl and it consists, at the second resonance, of the 
two TT-pulses producing swaps to the single-photon states, 
K2T4 — K,2Ts — hn, as well as the two 7r-pulses producing 
swaps to the two-photon states, \/2'«272 = \/2k2Tq = hn. 
The first and the third interference loops are included in 
the pulse sequence to annihilate the state superpositions, 
the corresponding phase shifts satisfying the relations. 



= TT, 05 — 6*7 = TT (mod27r). 



(54) 



E2 

h 



dt. 



The middle loop is required for generating a non-trivial 
gate. Specifically, the pulse sequence (T2 - Tg) produces 
on the states \geO) and \ggl) the following gate operation. 



f/(2)^vt/(^)55C/(-)53f/(^) = e--^^ J,, 

where Sn = diag[e~*^", 1]. For the levels l^el) and \gg2), 
the sequence produces the gate operation. 











Ti 



fj T3 T4 T5 Tg Tv 



t 



FIG. 6: Correct pulse sequence for performing a control- 
phase gate: time intervals Ti, T^, Tg and Tg are single-photon 
TT-pulses, whereas T2 and T% are two-photon yr-pulses; free 
evolution during times Tz,Tc, and Tt is added to annihilate 
excited photon states, and create non-trivial phase shift. 

Incorporating the modified gate operation into the se- 
quence, Eq. 152|) . we end up with the state evolution 
corresponding to a diagonal gate operation. 



I550) ^ 


I55O) - 


\ggO) ^ 


+ I.95O) 


\geO) ^ 


l^eO) - 


e-'^^\geQ) 




\egO) ^ 






|eeO) ^ 


\gei) - 


e-'«^|5el) ^ 


-e-*«5|eeO). 



The overall protocol produces the universal CPHASE 
gate (up to a common phase and single-qubit rotations), 



I55O) 
l5eO) 
leffO) 
|eeO) 



1 



1 



I5.9O) 

\geO) 
\egO) 

|eeO) 



(56) 
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For a specific choice of the angle, 65 = 7r/2(mod 27r), 
the gate in Eq. (|56l) can be further transformed to a 
CNOT gate using Hadamard rotations on qubit 2 as 
shown in Fig. [7| 



Qubit 1 



Qubit 2 



1 
i, 



H 



FIG. 7: Gate circuit for constructing a CNOT gate using tlie 
control-phase gate: a z-axis rotation is applied to qubit 1, and 
Hadamard gates H are applied to the second qubit. 



IV. DISCUSSION 

The protocol for the CPHASE gate considered in the 
previous section is much faster than the one suggested 
for a similar qubit-oscillator system and based on the 
dispersive qubit-oscillator coupling^- (by the ratio of the 
coupling frequency to the deviation from the exact res- 
onance (detuning)). On the other hand, the CPHASE 
gate protocol discussed in this paper is more complex 
and relatively slower than the protocol for direct qubit- 
qubit 2;2;-coupling considered in Refs. 12.33t the duration 
of the gate operation in the latter case is h/8 in the units 
of inverse coupling energy, while it is 2.7h for the proto- 
col presented in Fig. This illustrates the advantage of 
longitudinal, zz coupling (in the qubit eigenbasis), which 
is achieved for the charge qubits biased at the charge 
degeneracy point by current-current coupling^Si^ More 
common for charge qubits is the capacitive coupling, how- 
ever there the situation is different: this coupling has 
XX symmetry at the charge degeneracy point, and be- 
cause of inevitable difference in the qubit frequencies, 
the gate operation takes much longer time, prolonged by 
the ratio between the qubits frequency asymmetry and 
the coupling frequency. Recent suggestions to employ 
dynamic control methods to effectively bring the qubits 
into resonanca^^iS^ can speed up the gate operation. For 
these protocols, the gate duration is '--^ ft. in units of di- 
rect coupling energy, which is longer than in the case 
of zz coupling, but somewhat shorter than in our case. 
However, the protocol considered in this paper might be 
made faster by using pulse shaping. 

A specific feature of the present protocol is that there is 
always at least one qubit staying off-resonance: one qubit 
stays off-resonance during manipulations with the other 
qubit, and also the second qubit undergoes off-resonance 
excursions during the manipulation. During these peri- 
ods, free evolution of the qubits is assumed, however, a 
finite off-resonant qubit-cavity coupling violates this as- 
sumption and eventually negatively affects the fidelity 
of the gate. Therefore it is particularly important for 
this protocol to provide a weak off-resonance interaction. 



The latter is perturbatively estimated as ~ k'^/Ej, which 
is by factor, ~ Cc/Cs:, smaller than the resonant inter- 
action. In practice this factor would not exceed 1/100 
(coupling frequency ~ 80 MHz for qubit frequency > 5 
GHz), which would provide the fidelity of the gate not 
worse than one percent. 

To estimate a maximum number of qubits the circuit 
is able to accommodate, we assume reduced coupling 
strength, n/h ~ 10 MHz (which would still allow up to 
10 two-qubit operations during optimistic 1 /itsec deco- 
herence time). Assuming the same gate fidelity, 1%, and 
given the fact that the maximum qubit frequency is in 
the range of Ej/h ^ 10 GHz, we find that the num- 
ber of qubits can not be more than ten qubits. Thus 
we conclude that the resonant method of selective qubit 
addressing considered in the present paper allows to im- 
plement a small quantum processor suitable for testing 
the simplest quantum algorithms, although it has limited 
potential for larger scale applications. 
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APPENDIX A: CHARGE-PHASE REGIME 

In this appendix we derive Eq. H37() for the qubit-cavity 
interaction in the charge-phase regime. 

The starting point is the SCB Hamiltonian in Eq. H31|) 
at the degeneracy point, Ug = 1/2, written in the charge 
basis, \n). 



H = [Ec{n^l/2)'\n}{n\ 



n— — OQ 



- {Ej/2){\n + l){n\ + \n-l){n\)], (Al) 

We split the complete set of the charge eigenstates, 
—00 < n < 00, in the positive and negative charge sub- 
sets labelled with a ~1, 1, and m — ...,2,1 such that 



m 



n, n > 0, 
n < 0. 



m = 1 - n, n < 0. (A2) 
In the basis \m,a), the Hamiltonian H acquires the form. 



H 



A B 
B A 



where A is tridiagonal matrix. 



A 



Ec 
4 



—p 9 — p 
1 



2Ej 
Ec 



(A3) 



(A4) 



10 



and matrix B contains only a single element, 




1 



(A5) 



Similarly, the interaction Hamiltonian, corresponding to 

Eqs. (El, (EHIi, (Eni), 



Hint — 


f(<.+«')E( 

n—oc 


\n + 


l){n\- 


\n- 


l)(n|), (A6) 


takes the form in the (m. 


cr)-representation. 






Hint — (a 4 


at) 


C 

-iD - 


iD ' 

-C 




(A7) 


where 
















-i i 

~i 


' ^^2 





1 


• (A8) 



A Hadamard rotation, H, in cr-space. 



1 



(A9) 



takes the basis \m "f), \m J,) to 

|m±) = (1/\/2)(|to t) ± \m I)), and transforms the ma- 
trices in Eqs. ljA3|l and (|A7|) . to the form. 



The eigenstates of the SCB Hamiltonian, the Bloch 
states, which consist of superpositions of many charge 
states, can be found by independent unitary rotations in 
the a = + and a — — subspaces due to block-diagonal 
form of the matrix H in Eq. HA10|I . 

The qubits states are chosen as the two lowest energy 
Bloch states. The corresponding energy eigenvalues oc- 
cupy the bottom right corners of the diagonalized ma- 
trices, [W{A + B)U]ii and [V^A- B)V]ii, where 
and V{p) are appropriate unitary matrices. Indeed, this 
is obviously true for the charge regime {p <C 1), when 
these matrix elements correspond to superpositions of the 
n = and n = 1 charge states. When the Josephson en- 
ergy increases, which corresponds to increasing potential 
energy in the Bloch Hamiltonian (p ~ 1), the Bloch ener- 
gies change without crossing?^, and therefore the bottom 
right corner matrix elements remain to be the lowest en- 
ergy eigenvalues. 

In the qubit subspace, the interaction Hamiltonian in 
Eq. HA10|) takes the form. 



H 



-at) 



[W{C-iD)V]n 
[VHC + ID)U]^^ 

(All) 

Denoting the matrix element of the Hermitian matrix in 
Eq. I|A11|I as k/2 = /(p)k/2, where f{p) is a complex 
numerical coefficient, we arrive at the qubit-cavity inter- 
action on the form, 



Hi, 



(kct+ -I- K*a-) . 



(A12) 



H = A1 + Ba,, H„u^{a + a)){Ca^+Day). (AlO) 
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